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Physics prospects of UV-filtered overlap quarks 
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Some key features of the overlap operator with a UV-filtered Wilson kernel are discussed. The first part concerns 
spectral properties of the underlying shifted hermitean Wilson operator and the relation to the observed speedup 
of the overlap construction. Next, the localization of the filtered overlap and its axial- vector renormalization 
constant are discussed. Finally, results of an exploratory scaling study for mud,rns and f-K,fK are presented. 



1. INTRODUCTION 

On the conceptual level the quest for exact chi- 
ral symmetry at finite lattice spacing has been 
completed. The closely related domain-wall (in 
the limit L5 — > 00) ^ and overlap ^ approach 
and also the parametrized fixed point action (in 
the limit of perfect parametrization) ^ yield lat- 
tice fermions which satisfy the Ginsparg- Wilson 
relation with on-shell chiral symmetry jj] 

75^? -f i?75 = , 75=75(l-ii?). (1) 

p 

The salient features of such fermions include: 

• Additive mass renormalization is absent 0]. 

• There is an index theorem linking n+ — n_ 
(where n± is the number of zero-modes of 
the massless Dirac operator D with 75 — 
±1) to the gluonic topological charge q on 
fine enough lattices [S]- 

• The Nielsen-Ninomiya theorem is evaded Pl . 

• Using "chirally rotated" quark fields to con- 
struct e.g. pseudoscalar densities and axi- 
alvector currents (see below) the theory has 
only O(a^) cut-off effects 7 . 

On a practical level the question how one may 
implement any such scheme at bearable cost in 
terms of CPU time remains a topic of active re- 
search. Considering the massless overlap operator 



D 



p\^w-pDw-p) ^''^\ (2) 

p — ±j\isi~p the shifted Wilson operator 
(0 < p < 2) as an example, two avenues are being 



with D^^^^p — I?w 



pursued. One question is what is the numerically 
most efficient way to implement the (.)~^/^ oper- 
ator prescription in Q . The second strategy is to 
ask whether a slight modification of the formula- 
tion can reduce the computational burden. 



2. OVERLAP WITH UV-FILTERING 

An obvious idea is to do some "massage" to the 
Wilson kernel in the overlap prescription, since 
there is nothing specific to the (plain) Dw in @ ; 
any legal doubler-free fermion action will do fine. 

There is a long history of "designer actions" tai- 
lored to have a spectrum sufficiently close to the 
GW circle, such that the inverse squareroot in 
((2Jl would minimize the number of forward appli- 
cations of the kernel on a source vector [H] . These 
actions typically invoke tunable parameters and 
the improvement is achieved through additional 
couplings, i.e. a single row or column of D in- 
volves more than the 51 entries (for SU{3) and 
chiral repr.) of Dw- Hence the challenge is to as- 
sess the effect of fewer forward applications versus 
each application getting more expensive. 

In addition^, there is quite some experience 
with filtered actions (including the overlap) 
IIOIIIII2] where only the covariant derivative is 
replaced, 

Up{x)ij{x+fi)-tp{x) ^ Vf,{x)tp{x+ft)-ij{x) , (3) 

tantamount to a change by an irrelevant operator. 
Here, V^{x) is defined via an APE [TS|, HYP [TUI 



^Some of the "designer actions" do involve smeared links. 
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or stout-link [21 recipe^. The point is that this 
maintains the sparseness of the Wilson operator 
and yet a significant speedup can be achieved. 
The obvious concern is that the filtering might 
impair the localization properties of the final D. 

2.1. Spectral properties of |_ffvv, p| 

We now focus on how the speedup of the over- 
lap construction with UV-filtering can be under- 
stood in terms of the spectrum of the underlying 
hermitean Wilson operator ifw,-p = 75(^w — p)- 
We use the Wilson gauge action, and for technical 
details we refer to the original pubhcation [TK] . 

The numerical effort to construct the massless 
operator Q is in good approximation propor- 
tional to the condition number 1/e of |_ffw.-p|- 
Fig. 1 shows 1 /e on 25 quenched configurations 
at /3 = 6.0. Evidently, the condition numbers 
get dramatically reduced through any type of fil- 
tering. In practice, overlap calculations use the 
projection trick, i.e. the subspace pertinent to 
the lowest few eigenvalues is treated exactly |16| . 
The interesting news is that even after this trick 
has been applied, the condition number in the 
orthogonal complement is improved through fil- 
tering (right panel). The data are for p=l, with 
p>l the saving is slightly smaller. 

It has been argued that the spectral density of 
|i?vv,-p| getting too large will eventually make the 
overlap construction break down on coarse lat- 
tices ^7]. Fig. 2 shows the cumulative eigenvalue 
distribution (CED) of \Hw,-p\ with p=l on an 
extremely coarse {f3 — 5.66) and a fairly smooth 
(/3 = 6.0) lattice. The slope at the origin is the 
quantity of interest, p(0). The filtering clearly re- 
duces p{0), but by no means does it reduce it to 
zero. A point worth noticing is that the unfiltered 
CED on the coarse lattice has a rather different 
shape. At strong coupling the radii of the Wilson 
spectrum are considerably smaller than 1, and our 
choice to keep the shift parameter p=l fixed lets 
us "loose" the fermion. In summary, the overlap 
construction breaks down on too coarse lattices, 
but the good news is that (i) one would notice 
from the CED and (ii) the breakdown gets de- 
layed through filtering (/? = 5.66 seems fine if at 



The latter choice is mandatory in a dynamical HMC. 



least 1 filtering step is applied). 

Another question is what happens to the spec- 
tral density at weak coupling. Fig. 3 shows the 
CED of |i/w.-p| in 2D in a log-log representa- 
tion. With filtering a much higher fraction of the 
eigenvalues is near A = 1, and in this sense the 
picture is similar to what one would get if one 
were to use an approximately chiral kernel - in 
that case the "mobility edge" would be at A = 1 
and p(0) = 0. The high statistics study in 2D 
presented on the right suggests that p(0) decays 
exponentially in f3 and that filtering continues to 
reduce p(0) in the weak coupling regime. 

2.2. Localization properties of Dov 

It is evident from (|2Jl that Dov cannot be ul- 
tralocal, but the question is whether it is local, 
i.e. whether the couplings in Dov{x,y) would de- 
cay exponentially in |a;— (in some norm). This 
is important to guarantee that LQCD with Dq^ 
is in the right universality class. 

We measure the localization function 

/(di) = sup{||(i?ovry)(a:)||2 I ||a:-0||i=di} (4) 

with rj a normalized source vector at the origin. 
Here, di is the distance in the "taxi driver" met- 
ric. The localization v is the "effective mass" of 
/(di) for a di which is sufhciently far from the 
maximal one (2L in a L'^ box) to avoid finite vol- 
ume effects. Fig. 4 shows a scan of z/ as a function 
of the shift parameter p at /3 = 6.0. Without fil- 
tering the optimized p is near 1.4 ^S]. After just 
1 APE step the extremum shifts towards p ~ 1 
and with 1 HYP step it is near p ~ 0.8. This 
finding should not come as a surprise, since in 
the free theory a value fmax — 0.58 is realized 
for popt — 0.54. We consider this a strong ratio- 
nale for not tuning p to an "optimum" value at 
some standard coupling but for staying with the 
canonical choice p = 1. In passing we note that 
p= \ would also minimize the condition number 
of |iJw,-p| in the free field limit. 

Our choice to use the "taxi driver" distance was 
motivated by the standards in the literature [T^ . 
Fig. 5 shows f{di) in the free theory for several 
directions (on-axis and space-diagonals) together 
with the supremum. One notices that the slope 
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Figure 1. Condition number 1/e of |i7w,-p| at p=l on 25 quenched /3 = 6.0, 16 configurations without 
(left) and with projection of the 14 lowest eigenmodes (right). Modest filtering significantly reduces 1/e. 
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Figure 2. Cumulative eigenvalue distribution of |i/w,-p| at p=l for (3 — 5.66 (left) and f3 — 6.00 (right) 
with 0,1,3 steps of APE or HYP filtering. This reduces the slope and thus p(0), but the latter quantity 
remains non-zero. Evidently, the unfiltered operator at /3 = 5.66 is in the wrong universality class. 




none 
1APE 
3APE 



10- 
10"' 

10-' 

10"' 



10-^ lO-' 10" p 

Figure 3. Cumulative eigenvalue distribution (log-log plot) of |iJw,-p| in the Schwinger model {Nf — 0, 
16^, 13 — 3.2, p—l) with 0,1,3 filtering steps (left) and the pertinent p(0) as a function of (3 (right). 
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Figure 4. Localization v versus p at /3 = 6.0 without filtering and after 1 step of APE/HYP smoothing 
(left) and the same relation in the free theory (right); here the optimum p is significantly smaller than 1. 
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Figure 5. Fall-off pattern (vs. "taxi driver" distance d\) of the free overlap operator on-axis (open circles) 
and for several directions (other symbols) with supremum (full line) at p = 1 (left) and p = 0.54 (right). 
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Figure 6. AWI mass versus bare mass at /3 = 6.0 and p—1 for several filtering recipes (left). The quadratic 
fits are unconstrained and still go through zero, indicating good chiral symmetry. The inverse slope Za 
versus (3 (right) shows a much milder dependence with filtering than without, even with p^l kept fixed. 
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in the 2D-diagonal direction is roughly a factor 
•\/2 smaher than on-axis, in the SD-diagonal di- 
rection (the one which dominates the supremum) 
it is about a factor -v/S- In other words, if we had 
chosen the "beehne" metric ||.||2, then the maxi- 
mum locahzation in the free-field case would have 
turned out to be roug hly 0.58^3- 1. 

2.3. Axialvector renormalization for Dov 

In phenomenological studies theoretical uncer- 
tainties of latticc-to-continuum renormalization 
factors often limit the precision of the final an- 
swer. In this respect it seems important that such 
factors would be close to 1 or (at least) show a 
mild dependence on the gauge coupling. 

We determine Za via the axial Ward identity 
(AWI) . Specifically, we use the "chirally rotated" 
quark fields jZj to construct the pseudoscalar den- 
sity and the (naive) axialvector current 

P{X) = ^^(x)75[(l-;^Z?ov)^2](x) 



(5) 

(6) 



2p 

Ay,{x) = i/>i(2;)7^75[(l - — i:iov)V'2](a;) 

where below the fields (flavors) V'1.2 will be taken 
as solutions to the massive overlap operator 



/ am\ „ 



(7) 



with mass mi_2, respectively. The "rotating" op- 
erator [7] in H5I6|I is still the massless one. With 
these at hand one forms the ratio 



p(i, mi, 777,2) 



E^(P(x,i)P^(O,0)) 



(8) 



Ex(V4A4(x,i)P-(0,0)) 
which yields the sum of the AWI quark masses, 

p(t,m,,m2)*^mf^' + mf^' + 0{a'). (9) 



Fig. 6 shows the plateau value ^(777,1 -1-7712) for a 
variety of (7771,7772) combinations versus the sum 
mi+m2. Since each mi+m2 is realized in various 
combinations and we see no spread, we conclude 
that isospin is a perfect symmetry (in this ratio 
and within our statistical precision). Using 

p — const + 75— (Til -I- 7772) + const (7^1+7712)^ 



Za 



as our fit ansatz, we determine the offset at the 
origin and Za- The former is consistent with zero, 



and this means that we have indeed good chiral 
symmetry. How the Za factors would depend on 
the gauge coupling is shown on the right. Our 
data are for p = 1 and stem from |TKri9) . The 
data for p=lA come from j^Ol- This compilation 
shows that already a single HYP step at fixed 
p=l manages to tame the /3-dependence of Za in 
a more efficient way than carefully tuning p could 
possibly do. In fact, the Za values of the filtered 
overlap are so close to 1 that even a perturbative 
evaluation at the 1-loop level might work fine (as 
was done, for a different operator, in ,21 ). 

2.4. Technical summary 

Given the discussion in the two previous sub- 
sections it is clear that the UV-filtering does not 
only serve the purpose of speeding up the over- 
lap construction, it really improves the physics 
properties of the resulting Dqv Let us take the 
opportunity to highlight some technical points: 

(a) Filtering yields an 0(a^ )-redefinition of the 
overlap (at fixed p, asmcar, Mtor)- 

(b) The implementation effort is minimal; just 
evaluate _Dov on a smeared copy of the orig- 
inal gauge configuration. 

(c) The idea is well-suited for dynamical simu- 
lations with HMC 122], provided the stout- 
link recipe |14| is used. 

(d) The advantages are independent of how 
the (.)~^^^ prescription is implemented and 
naturally carry over to the domain-wall ^ 
or Moebius varieties. 

Regarding point (a) it should be clear that there is 
no conflict with the statement in '5'6' that "the" 
overlap yields a sound deflnition of "the" topolo- 
gical charge of a configuration U. It is true that 
for some configurations in a given ensemble the 
filtered overlap produces a mode excess number 

HYP-ltr(75a/?«7^[C/]) (10) 



q 



2p 



which differs from the unfiltered version 



1 , 

— tr 75 aD, 
2p 



none 
ov 



(11) 



It is important to notice that the same fact holds 
w.r.t. two unfiltered overlap varieties with un- 
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Figure 7. Mp versus mi + 7712 on our coarsest 
lattice (1 HYP operator) and the pertinent scaling 
plot for 2mud'ro and {rnud+TTis)ro (all filterings). 
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Figure 8. fp versus M|, on our coarsest lattice 
(1 HYP) and the scaling plot for /^ro and fx^o- 



equal p. The statement in |5I6) is that the num- 
ber of such "ambiguous" configurations dies out 
quickly enough with /3 ^ cx) to make the impact 
on any physical observable (e.g. Xtop) an O(a^) 
effect. These issues are discussed in jl2l24ll5j . 

3. EXPLORATORY SCALING TEST 

Given the promising localization properties of 
the filtered overlap, it is natural to ask whether 
this will lead to a better scaling behavior. There 
have been other scaling studies with the over- 
lap action [20124825] . but here we shall focus on 
the first continuum result for m„d, and /tt, 
from quenched overlap data |26j . 

On the technical level, one more ingredient is 
needed, Zs = Z:^ . We compute it with the RI- 



MOM method j27j. It turns out that the filtering 
leads to a rather mild /3-dependence of the scalar 
renormalization factor, too j26|. 

For the scaling study we use 4 couplings, (3 ~ 
5.66, 5.76, 5.84, 6.0, with lattice dimensions such 
as to have one matched spatial box size L = 1.5 fm. 
We choose 4 quark masses, ranging from about 
imPi^y" to mPl^y^ and 4 filterings (1 APE, 3 APE, 
1 HYP, 3 HYP). We do, however, not consider the 
unfiltered operator, simply because it is too de- 
manding for our computational resources. 

Fig. 7 shows the pseudoscalar mass squared ver- 
sus the sum of the valence quark masses (the Som- 
mer radius makes them dimensionless) on our 
coarsest lattice with the 1 HYP version of -Dov 
The analogous plots at the other couplings look 
rather similar (not shown). Imposing the phys- 
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ical value M^tq = 1.251 yields a definition of 
{rris + TUudji'o at this particular lattice spacing. 
Likewise, imposing A/^rg = 0.3537 yields 2raudTo 
at this coupling. Repeating this procedure for the 
other operators and on the finer lattices, we are 
finally in a position to extrapolate to the contin- 
uum. We obtain the quenched continuum values 

mf^(2GeV) = 119(10)(7) MeV (12) 

rris/mud = 23.7(7.1)(4.5) (13) 

with the first error being statistical, the second 
systematic (up to quenching). 

Fig. 8 shows the pseudoscalar decay constant 
versus the pseudoscalar mass squared, again on 
our coarsest lattice with the 1 HYP operator. Su- 
perimposing the data obtained with other filter- 
ings or at weaker coupling would reveal a rather 
good agreement. Using the same values for Mxro 
and M^rp as before, we get fK^Q and fT^ro at this 
coupling. Repeating this procedure on the finer 
lattices, we can extrapolate the pseudoscalar de- 
cay constants to the continuum. We obtain 

fK = 170(10)(2)MeV (14) 

Ik/U = 1.17(4)(2) (15) 

in the continuum with a systematic uncertainty 
that does not include quenching effects. 

Admittedly, our continuum extrapolations are 
somewhat courageous - we start from a lattice as 
coarse as = 1 GeV and have = 2 GeV on 
the finest one. Given this range and our statistical 
precision, we can only conjecture that we are in 
the Symanzik scaling regime. But in a strict sense 
this statement remains true with any statistical 
precision and any range of couplings. 

It is interesting to compare our continuum re- 
sults to those of • The central values are con- 
sistent, but the relevant aspect is the error. The 
CP-PACS result was thought to be the "final" 
quenched value and consumed about half a year 
of runtime on their then new machine. Our calcu- 
lation used an equivalent of less than 1 year on a 
stand-alone PC. And yet the net combined error 
(statistical and systematics, without quenching) 
is almost the same. This illustrates that the ab- 
sence of additive mass renormalization, due to ex- 
act GW symmetry, has important practical con- 
sequences. 



4. SUMMARY 

The main features of overlap fermions with UV- 
filtering may be summarized as follows. 

1. The overlap construction per se (with any 
undoubled kernel) brings a variety of invalu- 
able conceptual features like on-shell chi- 
ral symmetry and a sound definition of the 
topological charge via the index theorem. 

2. Using an UV- filtered Wilson kernel several 
technical advantages can be achieved: there 
is no need to tune the shift parameter p, the 
locality of the resulting Dov and the nor- 
mality of the kernel D^f^^-p are improved, 
furthermore H^^-p has a better condition 
number, and renormalization factors like 
Za depend only weakly on the coupling. 

3. The details of the filtering procedure seem 
rather irrelevant; one may use APE, HYP 
or stout-smearing with any set of parame- 
ters. The crucial issue from a conceptual 
viewpoint is that these parameters (cksmcar 
and Alitor) are kept fixed when (3 changes, 
just like the other parameter p must stay 
fixed. In practice it seems advisable to stay 
with a moderate amount of filtering, e.g. 1 
to 3 iterations with standard parameters. 

4. Our exploratory scaling study suggests that 
such filtered overlap quarks enjoy good 
scaling properties already on rather coarse 
(a^^ ~ 1 — 1.5 GeV) lattices. Under the pro- 
viso that more detailed studies with higher 
statistics confirm this finding, it seems con- 
ceivable that overlap quarks become the 
method of choice whenever high-precision 
results in the continuum are sought-for. 
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